{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}C:\Users\ondre\OneDrive\Plocha\RnP\Revisions\Dataverse\svr_jeps_replication_log3.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res}20 Jan 2023, 13:20:51

{com}. **************************************************************

. 
. *** Replication code for "The 'Commitment trap' Revisited: ***

. 
. *** Experimental Evidence on Ambiguous Nuclear Threats" by ***

. 
. *** Michal Smetana, Marek Vranka, and Ondrej Rosendorf     ***

. 
. **************************************************************

. 
. 
. 
. *** The code was written in Stata 17.0 BE-Basic Edition ***

. 
. 
. 
. *** Please reach out to ondrej.rosendorf@fsv.cuni.cz if you have any questions concerning this replication file ***

. 
. 
. 
. *** IMPORTANT: This file is accompanied by the svr_jeps_replication_data3 dataset ***

. 
. 
. 
. *** Before proceeding with the replication, please make sure that the "asdoc", "coefplot", "estout" and "betterbar" package is installed ***

. 
. 
. 
. *** To install the estout package, use the following command ***

. 
. 
. 
. ssc install asdoc, replace
{txt}checking {hilite:asdoc} consistency and verifying not already installed...
all files already exist and are up to date.

{com}. 
. 
. 
. *** To install the coefplot package, use the following command ***

. 
. 
. 
. ssc install coefplot, replace
{txt}checking {hilite:coefplot} consistency and verifying not already installed...
all files already exist and are up to date.

{com}. 
. 
. 
. *** To install the estout package, use the following command ***

. 
. 
. 
. ssc install estout, replace
{txt}checking {hilite:estout} consistency and verifying not already installed...
all files already exist and are up to date.

{com}. 
. 
. 
. *** To install the betterbar package, use the following command ***

. 
. 
. 
. ssc install betterbar, replace
{txt}checking {hilite:betterbar} consistency and verifying not already installed...
all files already exist and are up to date.

{com}. 
. 
. 
. *** Setting the output scheme to black and white ***

. 
. 
. 
. set scheme s1mono

. 
. 
. 
. ***************************************************

. 
. *** Replication of the results in the main text ***

. 
. ***************************************************

. 
. 
. 
. *** Figure 7 (bar graph) - Response likelihood across treatments ***

. 
. 
. 
. * Generating a labelled version of the scenario_n variable

. 
. recode scenario_n (0=0 "Ambiguous threat") (1=1 "Explicit threat"), generate(scenario_n_label)
{txt}(0 differences between {bf:scenario_n} and {bf:scenario_n_label})

{com}. 
. 
. 
. * Generating a bar graph with error bars

. 
. betterbarci economic conventional invade cyber nuclear nothing, over(scenario_n_label) bar format(%2.0f) vertical ylabel(0(10)100)

. 
. 
. 
. * Exporting the bar graph (Figure 7)

. 
. graph export F07.png
{txt}{p 0 4 2}
file {bf}
F07.png{rm}
saved as
PNG
format
{p_end}

{com}. 
. 
. 
. *************************************************

. 
. *** Replication of the results in Appendix 12 ***

. 
. *************************************************

. 
. 
. 
. *** Appendix 12, Table 1 - Independent samples t-test (Welch) ***

. 
. 
. 
. * Running the t-test (Welch) for individual likelihood estimates by scenario_n and generating tables

. 
. asdoc ttest nuclear, by(scenario_n) welch replace title(Table 1: Welch's t-test results) save(A12T01)

{txt}Two-sample t test with unequal variances
{hline 9}{c TT}{hline 68}
   Group{col 10}{c |}{col 16}Obs{col 27}Mean{col 35}Std. err.{col 47}Std. dev.{col 59}[95% conf. interval]
{hline 9}{c +}{hline 68}
       0 {c |}{res}{col 12}     75{col 22}    50.72{col 34} 3.624817{col 46} 31.39184{col 58} 43.49739{col 70} 57.94261
       {txt}1 {c |}{res}{col 12}     76{col 22} 54.42105{col 34} 4.047146{col 46}  35.2822{col 58} 46.35872{col 70} 62.48338
{txt}{hline 9}{c +}{hline 68}
Combined {c |}{res}{col 12}    151{col 22} 52.58278{col 34} 2.713735{col 46} 33.34694{col 58}  47.2207{col 70} 57.94487
{txt}{hline 9}{c +}{hline 68}
    diff {c |}{res}{col 22}-3.701053{col 34} 5.433111{col 58}-14.43672{col 70} 7.034618
{txt}{hline 9}{c BT}{hline 68}
    diff = mean({res}0{txt}) - mean({res}1{txt})                                      t = {res} -0.6812
{txt}H0: diff = 0                             Welch's degrees of freedom = {res} 149.387

    {txt}Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = {res}0.2484         {txt}Pr(|T| > |t|) = {res}0.4968          {txt}Pr(T > t) = {res}0.7516

{txt}{p 20 20 2}Table 1: Welch's t-test results{p_end}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:nuclear by scenario_n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:    50.72}}}{space 1}{space 1}{ralign 9:{res:{sf:   54.421}}}{space 1}{space 1}{ralign 9:{res:{sf:   -3.701}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:St_Err}{space 1}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:nuclear by scenario_n}}{c |}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:    5.433}}}{space 1}{space 1}{ralign 9:{res:{sf:      -.7}}}{space 1}{space 1}{ralign 9:{res:{sf:     .497}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}{space 1}{ralign 9:St_Err}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 13:nuclear by~n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:    50.72}}}{space 1}{space 1}{ralign 9:{res:{sf:   54.421}}}{space 1}{space 1}{ralign 9:{res:{sf:   -3.701}}}{space 1}{space 1}{ralign 9:{res:{sf:    5.433}}}{space 1}

{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}
{space 0}{res:{lalign 13:nuclear by~n}}{c |}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:      -.7}}}{space 1}{space 1}{ralign 9:{res:{sf:     .497}}}{space 1}
- mata: asdoctable("`save'", "`r(colnames)'", "`r(rownames)'", "T", "`r(ttitle)'", `fs', `abb', "`append'", "`noheader'", "`rowappend'", "," ,`frs',"", "`dec'", "`cmd'", "`label'", "`tzok'", "`font'", "`fhc'", "`fhr'", "`notes'", "`store'")
= mata: asdoctable("A12T01.doc", "obs1 , obs2 , Mean1 , Mean2 , dif , St_Err , t_value , p_value", "nuclear by scenario_n: 0 1,", "T", "Table 1: Welch's t-test results", 20, 20, "", "", "", "," ,1400,"", "3", "", "", "", "Garamond", "", "", "", "")
{res}{txt}{p 0 4 2}
(file {bf}
A12T01.doc{rm}
not found)
{p_end}
- set trace off
Click to Open File:  {browse "A12T01.doc"}

{com}. 
. asdoc ttest conventional, by(scenario_n) welch rowappend save(A12T01)
{txt}(File A12T01.doc already exists, option {bf:append} was assumed)

Two-sample t test with unequal variances
{hline 9}{c TT}{hline 68}
   Group{col 10}{c |}{col 16}Obs{col 27}Mean{col 35}Std. err.{col 47}Std. dev.{col 59}[95% conf. interval]
{hline 9}{c +}{hline 68}
       0 {c |}{res}{col 12}     75{col 22} 72.13333{col 34} 2.602794{col 46} 22.54085{col 58} 66.94715{col 70} 77.31951
       {txt}1 {c |}{res}{col 12}     76{col 22} 66.14474{col 34} 3.785091{col 46} 32.99766{col 58} 58.60445{col 70} 73.68502
{txt}{hline 9}{c +}{hline 68}
Combined {c |}{res}{col 12}    151{col 22} 69.11921{col 34}  2.30764{col 46} 28.35676{col 58} 64.55953{col 70} 73.67888
{txt}{hline 9}{c +}{hline 68}
    diff {c |}{res}{col 22} 5.988596{col 34} 4.593631{col 58}-3.096691{col 70} 15.07388
{txt}{hline 9}{c BT}{hline 68}
    diff = mean({res}0{txt}) - mean({res}1{txt})                                      t = {res}  1.3037
{txt}H0: diff = 0                             Welch's degrees of freedom = {res} 134.185

    {txt}Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = {res}0.9027         {txt}Pr(|T| > |t|) = {res}0.1946          {txt}Pr(T > t) = {res}0.0973

{txt}{p 20 20 2}Two-sample t test with equal variances{p_end}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:conventional by scen~n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:  72.1335}}}{space 1}{space 1}{ralign 9:{res:{sf:  66.1445}}}{space 1}{space 1}{ralign 9:{res:{sf:   5.9885}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:St_Err}{space 1}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:conventional by scen~n}}{c |}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:   4.5935}}}{space 1}{space 1}{ralign 9:{res:{sf:      1.3}}}{space 1}{space 1}{ralign 9:{res:{sf:    .1945}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}{space 1}{ralign 9:St_Err}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 13:convention~n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:  72.1335}}}{space 1}{space 1}{ralign 9:{res:{sf:  66.1445}}}{space 1}{space 1}{ralign 9:{res:{sf:   5.9885}}}{space 1}{space 1}{ralign 9:{res:{sf:   4.5935}}}{space 1}

{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}
{space 0}{res:{lalign 13:convention~n}}{c |}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:      1.3}}}{space 1}{space 1}{ralign 9:{res:{sf:    .1945}}}{space 1}
- mata: asdoctable("`save'", "`r(colnames)'", "`r(rownames)'", "T", "`r(ttitle)'", `fs', `abb', "`append'", "`noheader'", "`rowappend'", "," ,`frs',"", "`dec'", "`cmd'", "`label'", "`tzok'", "`font'", "`fhc'", "`fhr'", "`notes'", "`store'")
= mata: asdoctable("A12T01.doc", "obs1 , obs2 , Mean1 , Mean2 , dif , St_Err , t_value , p_value", "conventional by scenario_n: 0 1,", "T", "Two-sample t test with equal variances", 20, 20, "append", "header not needed", "rowappend", "," ,1400,"", "3", "", "", "", "Garamond", "", "", "", "")
{res}{txt}- set trace off
Click to Open File:  {browse "A12T01.doc"}

{com}. 
. asdoc ttest cyber, by(scenario_n) welch rowappend save(A12T01)
{txt}(File A12T01.doc already exists, option {bf:append} was assumed)

Two-sample t test with unequal variances
{hline 9}{c TT}{hline 68}
   Group{col 10}{c |}{col 16}Obs{col 27}Mean{col 35}Std. err.{col 47}Std. dev.{col 59}[95% conf. interval]
{hline 9}{c +}{hline 68}
       0 {c |}{res}{col 12}     75{col 22} 66.50667{col 34}  2.80369{col 46} 24.28066{col 58} 60.92019{col 70} 72.09314
       {txt}1 {c |}{res}{col 12}     76{col 22} 55.82895{col 34}  3.59358{col 46} 31.32811{col 58} 48.67017{col 70} 62.98773
{txt}{hline 9}{c +}{hline 68}
Combined {c |}{res}{col 12}    151{col 22} 61.13245{col 34} 2.316458{col 46} 28.46511{col 58} 56.55535{col 70} 65.70955
{txt}{hline 9}{c +}{hline 68}
    diff {c |}{res}{col 22} 10.67772{col 34} 4.557905{col 58} 1.668081{col 70} 19.68736
{txt}{hline 9}{c BT}{hline 68}
    diff = mean({res}0{txt}) - mean({res}1{txt})                                      t = {res}  2.3427
{txt}H0: diff = 0                             Welch's degrees of freedom = {res} 142.882

    {txt}Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = {res}0.9897         {txt}Pr(|T| > |t|) = {res}0.0205          {txt}Pr(T > t) = {res}0.0103

{txt}{p 20 20 2}Two-sample t test with equal variances{p_end}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:cyber by scenario_n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:  66.5065}}}{space 1}{space 1}{ralign 9:{res:{sf:   55.829}}}{space 1}{space 1}{ralign 9:{res:{sf:  10.6775}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:St_Err}{space 1}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:cyber by scenario_n}}{c |}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:    4.558}}}{space 1}{space 1}{ralign 9:{res:{sf:     2.35}}}{space 1}{space 1}{ralign 9:{res:{sf:    .0205}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}{space 1}{ralign 9:St_Err}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 13:cyber by s~n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:  66.5065}}}{space 1}{space 1}{ralign 9:{res:{sf:   55.829}}}{space 1}{space 1}{ralign 9:{res:{sf:  10.6775}}}{space 1}{space 1}{ralign 9:{res:{sf:    4.558}}}{space 1}

{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}
{space 0}{res:{lalign 13:cyber by s~n}}{c |}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:     2.35}}}{space 1}{space 1}{ralign 9:{res:{sf:    .0205}}}{space 1}
- mata: asdoctable("`save'", "`r(colnames)'", "`r(rownames)'", "T", "`r(ttitle)'", `fs', `abb', "`append'", "`noheader'", "`rowappend'", "," ,`frs',"", "`dec'", "`cmd'", "`label'", "`tzok'", "`font'", "`fhc'", "`fhr'", "`notes'", "`store'")
= mata: asdoctable("A12T01.doc", "obs1 , obs2 , Mean1 , Mean2 , dif , St_Err , t_value , p_value", "cyber by scenario_n: 0 1,", "T", "Two-sample t test with equal variances", 20, 20, "append", "header not needed", "rowappend", "," ,1400,"", "3", "", "", "", "Garamond", "", "", "", "")
{res}{txt}- set trace off
Click to Open File:  {browse "A12T01.doc"}

{com}. 
. asdoc ttest invade, by(scenario_n) welch rowappend save(A12T01)
{txt}(File A12T01.doc already exists, option {bf:append} was assumed)

Two-sample t test with unequal variances
{hline 9}{c TT}{hline 68}
   Group{col 10}{c |}{col 16}Obs{col 27}Mean{col 35}Std. err.{col 47}Std. dev.{col 59}[95% conf. interval]
{hline 9}{c +}{hline 68}
       0 {c |}{res}{col 12}     75{col 22} 69.58667{col 34} 2.565193{col 46} 22.21522{col 58} 64.47541{col 70} 74.69792
       {txt}1 {c |}{res}{col 12}     76{col 22}    61.25{col 34} 3.601906{col 46} 31.40069{col 58} 54.07463{col 70} 68.42537
{txt}{hline 9}{c +}{hline 68}
Combined {c |}{res}{col 12}    151{col 22} 65.39073{col 34} 2.234527{col 46} 27.45833{col 58} 60.97551{col 70} 69.80594
{txt}{hline 9}{c +}{hline 68}
    diff {c |}{res}{col 22} 8.336667{col 34} 4.421984{col 58}-.4076433{col 70} 17.08098
{txt}{hline 9}{c BT}{hline 68}
    diff = mean({res}0{txt}) - mean({res}1{txt})                                      t = {res}  1.8853
{txt}H0: diff = 0                             Welch's degrees of freedom = {res} 136.753

    {txt}Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = {res}0.9692         {txt}Pr(|T| > |t|) = {res}0.0615          {txt}Pr(T > t) = {res}0.0308

{txt}{p 20 20 2}Two-sample t test with equal variances{p_end}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:invade by scenario_n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:  69.5865}}}{space 1}{space 1}{ralign 9:{res:{sf:    61.25}}}{space 1}{space 1}{ralign 9:{res:{sf:   8.3365}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:St_Err}{space 1}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:invade by scenario_n}}{c |}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:    4.422}}}{space 1}{space 1}{ralign 9:{res:{sf:      1.9}}}{space 1}{space 1}{ralign 9:{res:{sf:    .0615}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}{space 1}{ralign 9:St_Err}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 13:invade by ~n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:  69.5865}}}{space 1}{space 1}{ralign 9:{res:{sf:    61.25}}}{space 1}{space 1}{ralign 9:{res:{sf:   8.3365}}}{space 1}{space 1}{ralign 9:{res:{sf:    4.422}}}{space 1}

{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}
{space 0}{res:{lalign 13:invade by ~n}}{c |}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:      1.9}}}{space 1}{space 1}{ralign 9:{res:{sf:    .0615}}}{space 1}
- mata: asdoctable("`save'", "`r(colnames)'", "`r(rownames)'", "T", "`r(ttitle)'", `fs', `abb', "`append'", "`noheader'", "`rowappend'", "," ,`frs',"", "`dec'", "`cmd'", "`label'", "`tzok'", "`font'", "`fhc'", "`fhr'", "`notes'", "`store'")
= mata: asdoctable("A12T01.doc", "obs1 , obs2 , Mean1 , Mean2 , dif , St_Err , t_value , p_value", "invade by scenario_n: 0 1,", "T", "Two-sample t test with equal variances", 20, 20, "append", "header not needed", "rowappend", "," ,1400,"", "3", "", "", "", "Garamond", "", "", "", "")
{res}{txt}- set trace off
Click to Open File:  {browse "A12T01.doc"}

{com}. 
. asdoc ttest economic, by(scenario_n) welch rowappend save(A12T01)
{txt}(File A12T01.doc already exists, option {bf:append} was assumed)

Two-sample t test with unequal variances
{hline 9}{c TT}{hline 68}
   Group{col 10}{c |}{col 16}Obs{col 27}Mean{col 35}Std. err.{col 47}Std. dev.{col 59}[95% conf. interval]
{hline 9}{c +}{hline 68}
       0 {c |}{res}{col 12}     75{col 22} 80.66667{col 34}  2.66903{col 46} 23.11448{col 58} 75.34851{col 70} 85.98483
       {txt}1 {c |}{res}{col 12}     75{col 22}       74{col 34} 3.282111{col 46} 28.42391{col 58} 67.46025{col 70} 80.53975
{txt}{hline 9}{c +}{hline 68}
Combined {c |}{res}{col 12}    150{col 22} 77.33333{col 34} 2.125685{col 46} 26.03422{col 58} 73.13295{col 70} 81.53372
{txt}{hline 9}{c +}{hline 68}
    diff {c |}{res}{col 22} 6.666667{col 34} 4.230364{col 58}-1.694997{col 70} 15.02833
{txt}{hline 9}{c BT}{hline 68}
    diff = mean({res}0{txt}) - mean({res}1{txt})                                      t = {res}  1.5759
{txt}H0: diff = 0                             Welch's degrees of freedom = {res} 143.934

    {txt}Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = {res}0.9414         {txt}Pr(|T| > |t|) = {res}0.1172          {txt}Pr(T > t) = {res}0.0586

{txt}{p 20 20 2}Two-sample t test with equal variances{p_end}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:economic by scenario_n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:  80.6665}}}{space 1}{space 1}{ralign 9:{res:{sf:       74}}}{space 1}{space 1}{ralign 9:{res:{sf:   6.6665}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:St_Err}{space 1}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:economic by scenario_n}}{c |}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:   4.2305}}}{space 1}{space 1}{ralign 9:{res:{sf:      1.6}}}{space 1}{space 1}{ralign 9:{res:{sf:     .117}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}{space 1}{ralign 9:St_Err}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 13:economic b~n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:  80.6665}}}{space 1}{space 1}{ralign 9:{res:{sf:       74}}}{space 1}{space 1}{ralign 9:{res:{sf:   6.6665}}}{space 1}{space 1}{ralign 9:{res:{sf:   4.2305}}}{space 1}

{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}
{space 0}{res:{lalign 13:economic b~n}}{c |}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:      1.6}}}{space 1}{space 1}{ralign 9:{res:{sf:     .117}}}{space 1}
- mata: asdoctable("`save'", "`r(colnames)'", "`r(rownames)'", "T", "`r(ttitle)'", `fs', `abb', "`append'", "`noheader'", "`rowappend'", "," ,`frs',"", "`dec'", "`cmd'", "`label'", "`tzok'", "`font'", "`fhc'", "`fhr'", "`notes'", "`store'")
= mata: asdoctable("A12T01.doc", "obs1 , obs2 , Mean1 , Mean2 , dif , St_Err , t_value , p_value", "economic by scenario_n: 0 1,", "T", "Two-sample t test with equal variances", 20, 20, "append", "header not needed", "rowappend", "," ,1400,"", "3", "", "", "", "Garamond", "", "", "", "")
{res}{txt}- set trace off
Click to Open File:  {browse "A12T01.doc"}

{com}. 
. asdoc ttest nothing, by(scenario_n) welch rowappend save(A12T01)
{txt}(File A12T01.doc already exists, option {bf:append} was assumed)

Two-sample t test with unequal variances
{hline 9}{c TT}{hline 68}
   Group{col 10}{c |}{col 16}Obs{col 27}Mean{col 35}Std. err.{col 47}Std. dev.{col 59}[95% conf. interval]
{hline 9}{c +}{hline 68}
       0 {c |}{res}{col 12}     75{col 22} 20.50667{col 34} 2.863273{col 46} 24.79667{col 58} 14.80147{col 70} 26.21186
       {txt}1 {c |}{res}{col 12}     76{col 22} 19.43421{col 34}  2.90392{col 46} 25.31578{col 58} 13.64931{col 70} 25.21912
{txt}{hline 9}{c +}{hline 68}
Combined {c |}{res}{col 12}    151{col 22} 19.96689{col 34} 2.032961{col 46} 24.98144{col 58} 15.94995{col 70} 23.98383
{txt}{hline 9}{c +}{hline 68}
    diff {c |}{res}{col 22} 1.072456{col 34} 4.078122{col 58}-6.985097{col 70}  9.13001
{txt}{hline 9}{c BT}{hline 68}
    diff = mean({res}0{txt}) - mean({res}1{txt})                                      t = {res}  0.2630
{txt}H0: diff = 0                             Welch's degrees of freedom = {res} 150.991

    {txt}Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = {res}0.6035         {txt}Pr(|T| > |t|) = {res}0.7929          {txt}Pr(T > t) = {res}0.3965

{txt}{p 20 20 2}Two-sample t test with equal variances{p_end}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:nothing by scenario_n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:  20.5065}}}{space 1}{space 1}{ralign 9:{res:{sf:   19.434}}}{space 1}{space 1}{ralign 9:{res:{sf:   1.0725}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}

{space 0}{hline 23}{c  TT}{hline 11}{hline 11}{hline 11}
{space 0}{space 0}{ralign 22:}{space 1}{c |}{space 1}{ralign 9:St_Err}{space 1}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 23}{c   +}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 23:nothing by scenario_n}}{c |}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 22:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:    4.078}}}{space 1}{space 1}{ralign 9:{res:{sf:      .25}}}{space 1}{space 1}{ralign 9:{res:{sf:     .793}}}{space 1}
{space 0}{hline 23}{c  BT}{hline 11}{hline 11}{hline 11}
{res}
{txt}{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:obs1}{space 1}{space 1}{ralign 9:obs2}{space 1}{space 1}{ralign 9:Mean1}{space 1}{space 1}{ralign 9:Mean2}{space 1}{space 1}{ralign 9:dif}{space 1}{space 1}{ralign 9:St_Err}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}{hline 11}
{space 0}{res:{lalign 13:nothing by~n}}{c |}{space 11}{space 11}{space 11}{space 11}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:       75}}}{space 1}{space 1}{ralign 9:{res:{sf:       76}}}{space 1}{space 1}{ralign 9:{res:{sf:  20.5065}}}{space 1}{space 1}{ralign 9:{res:{sf:   19.434}}}{space 1}{space 1}{ralign 9:{res:{sf:   1.0725}}}{space 1}{space 1}{ralign 9:{res:{sf:    4.078}}}{space 1}

{space 0}{space 0}{ralign 12:}{space 1}{c |}{space 1}{ralign 9:t_value}{space 1}{space 1}{ralign 9:p_value}{space 1}
{space 0}{hline 13}{c   +}{hline 11}{hline 11}
{space 0}{res:{lalign 13:nothing by~n}}{c |}{space 11}{space 11}
{space 0}{space 0}{ralign 12:0 1}{space 1}{c |}{space 1}{ralign 9:{res:{sf:      .25}}}{space 1}{space 1}{ralign 9:{res:{sf:     .793}}}{space 1}
- mata: asdoctable("`save'", "`r(colnames)'", "`r(rownames)'", "T", "`r(ttitle)'", `fs', `abb', "`append'", "`noheader'", "`rowappend'", "," ,`frs',"", "`dec'", "`cmd'", "`label'", "`tzok'", "`font'", "`fhc'", "`fhr'", "`notes'", "`store'")
= mata: asdoctable("A12T01.doc", "obs1 , obs2 , Mean1 , Mean2 , dif , St_Err , t_value , p_value", "nothing by scenario_n: 0 1,", "T", "Two-sample t test with equal variances", 20, 20, "append", "header not needed", "rowappend", "," ,1400,"", "3", "", "", "", "Garamond", "", "", "", "")
{res}{txt}- set trace off
Click to Open File:  {browse "A12T01.doc"}

{com}. 
. 
. 
. *** Appendix 12, Figure 1 (coefficient plot) - Likelihood of nuclear response (DV), no subset ***

. 
. 
. 
. * Running the OLS regression model (Model 1)

. 
. regres nuclear i.scenario_n

{txt}      Source {c |}       SS           df       MS      Number of obs   ={res}       151
{txt}{hline 13}{c +}{hline 34}   F(1, 149)       = {res}     0.46
{txt}       Model {c |} {res} 517.068916         1  517.068916   {txt}Prob > F        ={res}    0.4971
{txt}    Residual {c |} {res} 166285.646       149  1116.01105   {txt}R-squared       ={res}    0.0031
{txt}{hline 13}{c +}{hline 34}   Adj R-squared   ={res}   -0.0036
{txt}       Total {c |} {res} 166802.715       150   1112.0181   {txt}Root MSE        =   {res} 33.407

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}     nuclear{col 14}{c |} Coefficient{col 26}  Std. err.{col 38}      t{col 46}   P>|t|{col 54}     [95% con{col 67}f. interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
1.scenario_n {c |}{col 14}{res}{space 2} 3.701053{col 26}{space 2} 5.437325{col 37}{space 1}    0.68{col 46}{space 3}0.497{col 54}{space 4}-7.043173{col 67}{space 3} 14.44528
{txt}{space 7}_cons {c |}{col 14}{res}{space 2}    50.72{col 26}{space 2} 3.857479{col 37}{space 1}   13.15{col 46}{space 3}0.000{col 54}{space 4} 43.09757{col 67}{space 3} 58.34243
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{com}. 
. 
. 
. * Storing the estimates (Model 1)

. 
. estimates store M1

. 
. 
. 
. * Running the OLS regression model with controls (Model 2)

. 
. regres nuclear i.scenario_n i.male age income i.education_bin i.party

{txt}      Source {c |}       SS           df       MS      Number of obs   ={res}       150
{txt}{hline 13}{c +}{hline 34}   F(7, 142)       = {res}     7.08
{txt}       Model {c |} {res} 42899.5558         7  6128.50797   {txt}Prob > F        ={res}    0.0000
{txt}    Residual {c |} {res} 122845.278       142  865.107588   {txt}R-squared       ={res}    0.2588
{txt}{hline 13}{c +}{hline 34}   Adj R-squared   ={res}    0.2223
{txt}       Total {c |} {res} 165744.833       149  1112.38143   {txt}Root MSE        =   {res} 29.413

{txt}{hline 14}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}      nuclear{col 15}{c |} Coefficient{col 27}  Std. err.{col 39}      t{col 47}   P>|t|{col 55}     [95% con{col 68}f. interval]
{hline 14}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 1}1.scenario_n {c |}{col 15}{res}{space 2} 5.287279{col 27}{space 2} 4.858915{col 38}{space 1}    1.09{col 47}{space 3}0.278{col 55}{space 4}-4.317878{col 68}{space 3} 14.89244
{txt}{space 7}1.male {c |}{col 15}{res}{space 2}-29.16189{col 27}{space 2} 4.959354{col 38}{space 1}   -5.88{col 47}{space 3}0.000{col 55}{space 4}-38.96559{col 68}{space 3}-19.35818
{txt}{space 10}age {c |}{col 15}{res}{space 2} -.255744{col 27}{space 2}  .190346{col 38}{space 1}   -1.34{col 47}{space 3}0.181{col 55}{space 4} -.632022{col 68}{space 3}  .120534
{txt}{space 7}income {c |}{col 15}{res}{space 2}-1.137388{col 27}{space 2} .7245488{col 38}{space 1}   -1.57{col 47}{space 3}0.119{col 55}{space 4}-2.569684{col 68}{space 3} .2949083
{txt}1.education~n {c |}{col 15}{res}{space 2}-1.926298{col 27}{space 2} 5.437138{col 38}{space 1}   -0.35{col 47}{space 3}0.724{col 55}{space 4}-12.67449{col 68}{space 3} 8.821897
{txt}{space 13} {c |}
{space 8}party {c |}
{space 11}1  {c |}{col 15}{res}{space 2} 6.432551{col 27}{space 2} 6.193647{col 38}{space 1}    1.04{col 47}{space 3}0.301{col 55}{space 4}-5.811117{col 68}{space 3} 18.67622
{txt}{space 11}2  {c |}{col 15}{res}{space 2} 6.041501{col 27}{space 2} 6.082401{col 38}{space 1}    0.99{col 47}{space 3}0.322{col 55}{space 4}-5.982257{col 68}{space 3} 18.06526
{txt}{space 13} {c |}
{space 8}_cons {c |}{col 15}{res}{space 2} 80.36429{col 27}{space 2} 10.76853{col 38}{space 1}    7.46{col 47}{space 3}0.000{col 55}{space 4} 59.07694{col 68}{space 3} 101.6516
{txt}{hline 14}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}
{com}. 
. 
. 
. * Storing the estimates (Model 2)

. 
. estimates store M2

. 
. 
. 
. * Generating the coefficient plot (Appendix 12, Figure 1)

. 
. coefplot M1, bylabel(Model 1) || M2, bylabel (Model 2) ||, drop(_cons) xline(0) coeflabels(1.scenario_n = "{c -(}bf:Treatment (explicit – ambiguity){c )-}" 1.male = "Gender (male)" age = "Age" income = "Income" 1.education_bin = "Education (university degree)" 1.party = "Party (Democrat – Republican)" 2.party = "Party (Independent – Republican)")
{res}
{com}. 
. 
. 
. * Exporting the coefficient plot (Appendix 12, Figure 1)

. 
. graph export A12F01.png
{txt}{p 0 4 2}
file {bf}
A12F01.png{rm}
saved as
PNG
format
{p_end}

{com}. 
. 
. 
. *** Appendix 12, Table 2 - OLS regression (Model 1 and 2) ***

. 
. 
. 
. * Generating a table with results (Appendix 12, Table 2)

. 
. esttab M1 M2 using A12T02.rtf, noeqlines eqlabels(none) eform nogaps se r2 varlabels(1.scenario_n "Treatment (explicit - ambiguity)" 1.male "Gender (male)" age "Age" income "Income" 1.education_bin "Education (university degree)" 1.party "Party (Democrat - Republican)" 2.party "Party (Independent - Republican)" _cons "Constant") drop(0.scenario_n 0.male 0.education_bin 0.party) mtitle("Nuclear response likelihood" "Nuclear response likelihood") title(Table 2: OLS regression of nuclear response likelihood) nonumbers mlabels("Model 1" "Model 2")
{res}{txt}(output written to {browse  `"A12T02.rtf"'})

{com}. 
. 
. 
. ************************************************

. 
. *** Continue with svr_jeps_replication_code4 ***

. 
. ************************************************

. log close
      {txt}name:  {res}<unnamed>
       {txt}log:  {res}C:\Users\ondre\OneDrive\Plocha\RnP\Revisions\Dataverse\svr_jeps_replication_log3.smcl
  {txt}log type:  {res}smcl
 {txt}closed on:  {res}20 Jan 2023, 13:21:17
{txt}{.-}
{smcl}
{txt}{sf}{ul off}